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Time-dependent Internal DFT formalism and Kohn-Sham scheme. 
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We generalize to the time-dependent case the stationary Internal DFT / Kohn-Sham formalism 
presented in Ref. [l4| . We prove that, in the time-dependent case, the internal properties of a 
self-bound system (as an atomic nuclei or a Helium droplet) are all defined by the internal one-body 
density and the initial state. We set-up a time-dependent Internal Kohn-Sham scheme as a practical 
way to compute the internal density. The main difference with the traditional DFT / Kohn-Sham 
formalism is the inclusion of the center-of-mass correlations in the functional. 
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I. INTRODUCTION. 

Traditional Density Functional Theory (DFT) 0, 0, Q 
and its time-dependent generalization [J, |5j have evolved 
into standard tools for the description of electronic prop- 
erties in condensed-matter physics and quantum chem- 
istry through the simple local density instead of the less 
tractable TV-body wave function. Stationary DFT is 
based on the Hohenberg-Kohn (HK) theorem [||, which 
proves that, for any non-degenerate system of N Fermions 
or Bosons [1| put into a local external potential, the N- 
body ground-state wave function can be written as a 
functional of the local ground-state density. A similar 
theorem exists for the time-dependent case [1, Q , where 
a dependence on the initial state appears. The Kohn- 
Sham (KS) scheme 0] and its time-dependent general- 
ization 0, Q provide a straightforward method to com- 
pute self-consistently the density in a quantum frame- 
work, defining the non-interacting system (i.e. the local 
single-particle potential) which reproduces the exact den- 
sity. 

Traditional DFT is particularly well suited to study 
the electronic properties in molecules Q. As a molecule 
is a self-bound system, the corresponding Hamiltonian 
is translationally invariant (which ensures Galilean in- 
variance of the wave function Q ) , and one can apply the 
Jacobi coordinates method. This permits to decouple the 
center-of-mass (cm.) properties from the internal ones, 
and to treat correctly the redundant coordinate problem 
(i.e. the fact that one coordinate is redundant for the 
description of the internal properties (Io| ) and the cm. 
correlations. But as the nuclei are much heavier than the 
electrons, we can apply the Jacobi coordinates method to 
the nuclei only, so that only the nuclei will carry the cm. 
correlations, and use the clamped nuclei approximation. 
Then, one recovers the "external" potential of traditional 
DFT, of the form YliLi v ext{?i), which accounts for the 
nuclear background as seen by the electrons in the frame 
attached to the cm. of the nuclei. Thus, traditional 
DFT is particularly adapted to the study of the electronic 
properties in molecules 0]. It is implicitely formulated 
in the nuclear cm. frame [Tl| and the energy functional 
does not contain any cm. correlations. Of course, con- 



trary to the whole molecule, the pure electronic system is 
not a self-bound system: the v ex t potential breaks trans- 
lational invariance and is compulsory in order to reach 
bound states in the stationary case. 

For other self-bound systems, as isolated atomic nuclei 
or He droplets, the situation is intrinsically different be- 
cause the masses of all the particles (Fermions or Bosons) 
are of the same order of magnitude. As a consequence, 
to decouple the cm. properties from the internal ones, 
one has to apply the Jacobi coordinates method to all 
the particules. The redundant coordinate problem (thus 
the cm. correlations) will now concern all the particles 
and should be treated properly. If a DFT exists, the 
cm. correlations should be taken into account in the 
functional. 

Moreover, no "external" potential of the form 
X)<=i v ext(?i) can be justified in the corresponding self- 
bound Hamiltonians (we denote r; the N particules coor- 
dinates related to any inertial frame as the laboratory). 
One may be tempted to formulate a DFT using the tra- 
ditional DFT conclusions in the limit v ext — > 0, but this 
would lead to false and incoherent results because: 

• in the stationary case, the Hohenberg-Kohn theo- 
rem is valid only for external potentials that lead 
to bound many-body states |l3[, which is not the 
case anymore at the limit v ext — > for translational 
invariant particle-particle interactions [Til ]: 

• the form of v ex t is not translationally invariant, 
but translational invariance is a key feature of self- 
bound systems HEBE!; 

• traditional DFT concepts as formulated so far are 
not applicable in terms of a well-defined internal 
density pint, i.e. the density relative to the s yste m's 
cm. , which is of experimental interest [H, Il4 , flH | 
(it is for example measured in nuclear scattering 
experiments). 

Instead of the traditional DFT potential Y^iLi v ext( r i), 
one might be tempted to introduce an arbitrary transla- 
tional invariant potential of the form J2<=i "m((ri — R), 
where R = YljLi r j * s the totaL cm. of the particles. 
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This potential is an "internal" potential, i.e. is seen in 
the cm. frame, and in [Tij we underlined that it is the 
only form which satisfies all the key formal properties. 
However, i>j„ t should be zero in the purely isolated self- 
bound case. This is why in [lij we presented it as a 
mathematical "auxiliary" to reach our goal and showed 
that it can be dropped properly at the end, conserving 
all the conclusions. Through it (and using the Jacobi 
coordinates), we proved, by a different way than those 
found in 0, El, the stationary "Internal DFT" theo- 
rem: the internal many-body state can be written as a 
functional of pint- Then we formulated rigorously the 
corresponding "Internal" KS scheme (in the cm. frame). 
The main interest of this work is to give a first step to- 
wards a fundamental justification to the use of internal 
density functionals for stationary mean-field like calcula- 
tions of nuclei [ll[ or He droplets with effective inter- 
actions, showing that there exists an ultimate functional 
which permits to reproduce the exact internal density, 
which was not clear up to now. 

It is to be noted that the stationary Internal DFT / 
KS formalism gives a more fundamental justification than 
the Hartree-Fock (HF) framework to the stationary nu- 
clear mean-field like calculations. Indeed, HF does not 
contain quantum correlations, nor treats correctly the 
redundant coordinate problem, which introduces a spu- 
rious coupling between the internal properties and the 
cm. motion [Tel [20j. A way to overcome this problem in 
the stationary case is to perform projected HF (projec- 
tion before variation on cm. momentum), which permits 
to restore Galilean invariance, but at the price of aban- 
doning the independent-particle description [13, [HI, UU ■ 
Within the Internal DFT / KS formalism, we proved that 
the cm. correlations can be included in the energy func- 
tional / the KS potential [l4| , so that there would be no 
need for a cm. projection if the ultimate functional was 
known. 

It is a question of interest to generalize the station- 
ary Internal DFT / KS formalism to the time-dependent 
case. It would provide a first step towards a fundamental 
justification to the use of density functionals in nuclear 
time-dependent calculations with an effective mean-field 
[HI, [HI' an d would prove that the cm. correlations can 
be included in the functional. This last point is even 
more interesting that the spurious cm. motion problem 
remains in time-dependent HF [2l|, [22| , but that then the 
projected HF method becomes unmanageable and is not 
used in practise [22T |. 

In this paper, we propose to set up the time-dependent 
Internal DFT / KS formalism. The paper is organized as 
follows: we first apply the Jacobi coordinates method to 
the time-dependent full many-body Hamiltonian to de- 
couple the internal properties from the cm. ones, and 
define some useful "internal" observables, including the 
internal density (section II); then we show that the in- 
ternal many-body wave function (and thus the "inter- 
nal" mean values of all the observables) can be written 
as a functional of the internal density (section III); fi- 



nally, we develop the associated time-dependent Internal 
KS scheme as a practical scheme to compute the internal 
density (section IV). 

II. TIME-DEPENDENT N BODY 
FORMULATION. 

A. General formulation. 

In the time-dependent domain, the introduction of an 
explicitely time-dependent internal potential of the form 

N 

^2 v m t(vi — R; t) (1) 

?:=i 

takes a true meaning. This is because self-bound systems 
are plagued by a cm. problem. For instance, in the sta- 
tionary case, the cm. will be dclocalized in the whole 
space for isolated self-bound systems [H, [2 [HI- This 
does not occur in experiments because experimentally 
observed self-bound systems are not isolated anymore 
(they interact with the piece of matter they are inserted 
in which localizes the cm.). In the time domain, the 
cm. motion remains uncomparable to the experimental 
one (this will be discussed in more detail later), so that 
it would not make sense to introduce a time-dependent 
potential which would act on the cm. motion. It are 
the internal properties which are of true experimental 
interest (experimentalists always deduce those proper- 
ties (1B|). This justifies the introduction of an explicitely 
time-dependent potential of the form (fTJ), which would 
act on the internal properties only, and models the in- 
ternal effect (only) of time-dependent potentials used in 
experiments. Such a potential does not appear any more 
simply as a mathematical auxiliary (as for the station- 
ary Internal DFT / KS) and should not necessarily be 
dropped at the end. 

We thus start from a general translationally invariant 
TV-body Hamiltonian composed of the usual kinetic en- 
ergy term, a translationally invariant two-body potential 
it, which describes the particle-particle interaction, and 
an arbitrary translationally invariant " internal" potential 
Vi n t which contains an explicit time dependence 

N 2 N N 

H = E £ + E u ( r > v i) + E - R; t) . (2) 

i—l i,j=i i—1 

i>j 

For the sake of simplicity we assume a 2-body interaction 
u and N identical Fcrmions or Bosons. The generaliza- 
tion to 3-body etc interactions is straightforward; the 
generalization to different types of particles is underway. 

We rewrite the Hamiltonian ^ using the (N — 1) Ja- 
cobi coordinates {£ a ;c* = 1, ■ . ■ , N — 1} and the cm. 
coordinate R, defined as 

fi =r 2 - ri,f 2 = r 3 — , 
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N 



ZN-l — 



1 N 

AT 



N - 1 

N 

r 



(rjv-R), 



N 



(3) 



The £ Q are relative to the cm. of the other 1, . . . , a — 1 
particles and are independent from R. They are to be 
distinguished from the N "laboratory coordinates" r», 
and the N "cm. frame coordinates" (r* — R) relative to 
the total cm. R. As the {r^ — and the {r^ — R} 

can be rewritten as functions of the £ Q (in Appendix [A] 
is given the expression of the {r^ — R} as a function of 
the {£ Q } coordinates), the interaction u and the internal 
potential Vi nt can be rewritten as functions of the £ Q . We 
denote U and V the interaction potential and the internal 
potential in the Jacobi coordinates representation: 



N 



i>j 

N 

y~]vint(*i - R;*) 



f/(a,-,^-i) 

(4) 



i=l 



Of course we have U[u] and V^[wi„t]. The V[vj n t] potential 
is (N — 1) body in the Jacobi coordinates representation 
and cannot be written in a simple form in this represen- 
tation (see Appendix [A]) . Moreover, various vi n t can lead 
to the same V, which we will develop later. 

After having defined the conjugate momenta of R and 
£ Q , we can separate ([2]) into H = Hcm + -Hint, where 
(AI = Nm is the total mass) 



h 2 A 



R 



CM 



2AI 



(5) 



is a one-body operator acting in R space only, and (r a 



is the conjugate momentum of £ a and p a 
corresponding reduced mass) 



J a+l 



the 



hit 



N-l 9 
= V — 



+ f/[u](ei,...,^-l) 

+ ^(ft £AT-i;t) (6) 



is a (N — 1) body operator in the space. It contains 
the interaction and the internal potential. 

In the time-dependent case, we can choose freely the 
initial state i/j( t i> ■ ■ ■ > r N> *o)- We start from an initial 
state which can be written 

V>(ri, . . .,rjv;*o) = T(R;t ) ipint(Zi, ■ ■ .,^iv-i;*o) (7) 

in the Jacobi coordinates representation. This form does 
not mix the cm. motion with the internal one (mixing 
them would not make sense because the cm. motion 
does anyway not correspond to the experimental one) and 
corresponds to the form of the stationary state [lO, 13 • 
As Hcm and Hi nt act in two separate subspaces, the R 



and spaces (which implies [Hcm, H int } — 0), it is 
easy to show that the state \i/j(t)) can be built at all time 
t > to as a direct product of the form 

^(ri,...,rjv;i) = r ( R ;*) V>mt(£i: • • • j &v-i;*), (8) 
with 



H CM \T(t)) = ihdt\T(t)) 

Hi n t\^int(t)) = ihd t \4>int(t)) 



(9) 
(10) 



Hence, the Af-body wave function ijj can be separated 
into a one-body wave function T that depends on the 
position R of the cm. only, and an "internal" (N — 1) 
body wave function i/i mi that depends on the remaining 
(A — 1) Jacobi coordinates £ a . Of course, ipi n t could also 
be written as a function of the N laboratory coordinates 
Ti, but one of them would be redundant. T is solution of 
the free Schrodinger equation and describes the motion 
of the isolated system as a whole in any chosen inertial 
frame of reference (as the laboratory) . If one starts from 
a normalizable initial state |r(i )), |r(t)) is condemned 
to spread more and more. In the stationary limit, the 
only solutions of Eq. @ are plane waves, which are in- 
finitely spread (thus not normalizable). This does not 
correspond to experimental situations, where the system 
is not isolated anymore: interactions with other systems 
of the experimental apparatus localize the cm. But the 
formal decoupling between the cm. motion and the in- 
ternal properties obtained when using the Jacobi coordi- 
nates method permits to let the cm. motion to the choice 
of experimental conditions, the internal properties being 
comparable to the experimental ones. 



B. Some useful definitions. 

We define some quantities and relations that will be 
useful for the next considerations. In [3, HE HH is de- 
fined the internal one-body density 

Pmt (v,t)/N (11) 
= J dr\ ■ --dr N 5(R)|^ irit (ri, . . . , r N ; t)\ 2 5(r - (n - R)) 

= ( N _ 1 ) J d ^ ■ ■ -<%n-2 |V>mt(6, ■ ■ ■ ,&v-2, i^t;*)| 2 - 

It is is normalized to N. The laboratory density p(r, t) is 
obtained by convolution of pi n t with the cm. wave func- 
tion (following [Ulli): p(r,t) = J dR|T(R, t)\ 2 p int (r - 
R,t). 

We also introduced in [1J] the local part of the two- 
body internal density matrix 



7i„t(r, r';t) 
= Jdr 1 ■ --dr N tf(R)|V>tnt(ri 5 • ■ ■ ,i"Ar;t)| 2 

x<J(r-(r;-R))<5(r'-(r j¥ ;-R)) 



(12) 
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N(N - 1) fN - 1 
N-2 



N 



N - 1 



I ( t C r' + (jV-l)r jVr' ,\ 
Vint Kb ■ ■ ■ ,?AT-3, ]v=2 ' W^T' 1 ) 



It has the required normalisation to N(N — l)/2. Fol- 
lowing similar steps than in [H, [2^| , we can show that 
the local part of the two-body laboratory density matrix 
7(r, r', t) is obtained by convolution of 7i„t with the cm. 
wave function: j(r,r']t) = j c?R|r(R, t)\ 2 ji n t(r — R, r' — 
R;i). 

The definitions of pi nt (r,t) and 7j„ t (r, r';i) show 
clearly that they are defined in the cm. frame, i.e. that 
the r, r' coordinates are measured in the cm. frame (see 
the 5 relations in Ijlip and (fT!?)) ). Compared to the tra- 
ditional definitions, a <5(R) appears in the definition of 
the internal densities calculated with ipi nt in {r^} coordi- 
nates. As one of them is redundant, the <5(R) represents 



the dependence of the redundant coordinate on the oth- 
ers Hj. 

Another quantity that will be very useful is the one- 
body internal probability current, defined in Appendix [Bl 
(c.c. denotes the complex conjugate) 

jint(r,t)/N (13) 

fl / N \ 3 f 
\ N ~ 1 j / d£ 1 ...d£ N -2il>t nt (€i, ...,£ N -2,v,t) 



2mi 



xV„V>ini(£l, ■■;£,N-2,V]t) 



+ c.c. 



which satisfies the "internal" continuity equation 

dtPint(r,t)+V.j in t(r,t) =0. (14) 
Using (fl"3"|) . ([6]) and (JTOj) , we obtain the relation 



i-Q^int{v,t) 

N ( N 



(15) 



iv / iv \ J / r 
^~{ N _ l ) / d £i" d 6v-2 jV^mtCCi, ...,£ N -2,v;i)ihd t ip int (€i, ...^N-2,v;t) 

+Vu(ihd t ll>int(£l, £N-2,V,t))lp* nt (£l, ...^ N -2,V\t) + C.C.j 
T7 7 / d£ 1 ...d£, N _ 2 \ V^Vmt^l, -,v;t)- -i>* in t{£l, -,V,t) ~ Vint (6 . ^5 *) V„- "V'mtCCl, ^5 *) 

+V>*nt(£i> ■■■,z';i)V 1 ,(y[u](£i, ...,!/) + V[ui n i](^i, ...,t';*))V'mt(6) ■■■,^;i) + C.C.| 



iV / AT 
2mi 



which will be a key equation for the next considerations. tentials V[vi n t](£,i, —,£n-i', t) defines two maps [1, [f| 



III. TIME-DEPENDENT INTERNAL DFT 
THEOREM. 

A. Preliminaries. 

To prove the time-dependent Internal DFT theorem, 
we adapt the considerations of [1, Q to the internal 
Schrodinger equation p0[) . The main differences lie in 
the definition of the corresponding internal density (jlip 
and probability current (|13p . and in the fact that the 
potential V[vi n t\{^\, ...,£jv-i;t) cannot be written as the 
sum of one-body potentials in the Jacobi coordinates rep- 
resentation (which introduces some subtleties due to the 
cm. correlations and will bring us to use the integral 
mean value theorem to reach our goal). 

In what follows, we consider a given type of Fcrmions 
or Bosons, i.e. a given particle-particle interaction u. 
Solving the "internal" Schrodinger equation (fTUj) for a 
fixed initial state \ipint(to)) and for various internal po- 



F :Vlv in t](£x,...,£ N -i;t) \y in t{t)) 

G : V[v int \(€i, ...,£iv-i;i) -> Pint(r,t). (16) 

We first notice that two potentials Vi n t and v' int 
which lead to two potentials V[vi nt ](£,i, ...,£,N~i',t) and 
V[v' int }(£,i, ...,^N-i',t) that differ by a scalar function of 
time only c(t) , will give two wave functions that differ by 
a phase e~ la ^/ h only @,[5|: 

V[v' mt ] (a , Cjv-i ; t) - V [v mt ] (a, 6r-i; t) = c(t) 
=> W int {t)) = e- ia ^ n \i> int {t)l 

with d(t) = c(t). (17) 

Then, \ipint(t)) and \i^' int (t)) will give the same density 
Pmt(r,t) = p' int (r,t). The consequence is that the map 
G is not fully invertiblc. 

Let us discuss a bit about the condition (TT7]) . The form 



flU for V[v in t] implies V[v' mtl 
We define 



V\v v 



Av int (r ; t) = v' mt (r; t) - v int (r; t). 



(18) 
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It is to be noted that the condition Aii, nt (r;i) ^ c(t)/N 
is necessary but not sufficient to ensure the condition 
(fT7|) . which can be rewritten V[Avi nt ](^i, CiV-i! t) ^ 
c(t). Indeed, it is possible to have A-y irat (r;t) ^ c(t)/N 
and nevertheless V[Av int }(^i, ...,^N-i',t) = c(t), because 
compensations due to the cm. correlations can happen. 

Let us reason on the two particulcs case, where only 
one Jacobi coordinate is sufficient to describe the internal 
properties. We have (see Appendix [K§: V[Avi nt ](£,i; t) = 
Av int (-^i;t) + A<u in t(|£i;i)( = Y?i=i v v&{*i -R-!*))- 
We see that if Avi n t(r; t) is an odd function of r at all t 
(up to an additional time-dependent function), we have 
V[Av int \(£i;t) = c(t) => p int = p' int - This is due to 
the cm. correlations, that the non-trivial form of V re- 
flects. If Avi n t tends to move the first particule in one 
direction, the second particule will tend to move in the 
opposite direction because of the cm. correlations. But 
if this potential counter-acts perfectly the motion of the 
second particule (as does an odd potential in the cm. 
frame), then the particules remain stuck and the density 
unchanged. 

The same can occur for an arbitrary number of partic- 
ules. For instance, as J2iLi( r i~^-) = 0, it is obvious with 
flU) and (HHI) that every Av int {r;t) = b(t).r + c(t)/N will 
yield U[At>i n t] = c(t) (even if this form for Avi n t leads to 
internal potentials which are not null at infinity). Again, 
this is because if a potential counter acts perfectly the 
motion due to the cm. correlations, the particules re- 
main stuck and the density unchanged. In what follows, 



we consider only internal potentials Vi nt and v' int that 
lead to V[Av lnt ] ^ c(t). 

We come back to Eq. (fT7|) and denote \ipi n t(t)) = 
e -lQ! ( t )/ fi |'0° nt (t)) where we define tp° nt as the wave func- 
tion obtained for the choice c(t) = 0, i.e. associated to a 
V[vi n t}(£,i, ...,^N-i',t) where no additive time-dependent 
function can be split. If we prove that the map G 
is invertible up to an additive time-dependent function 
c(i), then tpf nt is fixed by pi„t through the relation 
lCt(*)) = ^G-Vmt(r,t), which implies that \^ nt {t)) 
can be written as a functional of the internal density 
Pint defined in (jlip . Consequently, any expectation value 
of an operator O which does not contain a time deriva- 
tive can be written as a functional of the internal den- 
sity (as the phase cancels out): (V'mt (*)|0| V'int (*)) = 

(C*M(*)|oK*W(t)). 

We thus have to show that a propagation of 
(JTU]) with two potentials Vi n t and v' int that yield 
V[Avint]{£,i, ...,£jv-i;t) 7^ c {t) will produce two different 
internal densities pi nt and p' int . 



B. The proof. 

We start from a fixed initial state |Vw(io)) and prop- 
agate it with two with two potentials Vi nt and v' int that 
give V[Avi n t}(£i, ...,£jv_i;t) ^ c(i). We deduce from Eq. 



i^-(hnt(r,t) -j' mt (r,t)) =—(^—\ I d6v_ 2 |Vint(Ci, ■■■,v;t Q )\ 2 V y V[Av lnt ](i 1 ,...,v; t ) 
at \ J t=t mi\N — 1/ J 



.(19) 



Using the "internal" continuity relation (jT4j) we obtain 

(f>int(r,t) -Pintfatj) t=f = 

Nf N \ 3 f 

m \ N - l ) Vr ' / ^ 1 - c ^ JV - 2 l^ int ^ 1 '" 



(20) 



N 



N-l 



r;t )\ 2 V„V[Av m t}{Zi,...,v;t ) 



We now make the only hypothesis which is used in this derivation. Following [J, [f| we restrict the set of potentials 
Vint to those that can be expanded into Taylor series with respect to the time at the initial time to (which is a 
reasonable hypothesis for physical potentials). As we supposed that T^[Auj„t](£i, £jv-i; t) ^ c(t), we have (k is a 
positive integer) 



V[Av int ](€i, ...,^N-V,t) c(t) => 3k : Wk(£i, ...,£jv-i; *o) constant, 



(21) 



where 



wkdx, ...,^N-i;t ) = QjpV[Avi nt ](€i, ...,^N-i;t) 



It is to be noted that the condition ^-A!)j„f(r; t) ^ constant => V r JprAvi„ t (r; t) 

t — to 

the condition (|21[) (see (0| and ([22])), but not sufficient. 



t=t 



(22) 

7^ is necessary to ensure 
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In what follows, we consider k as the smallest positive integer such that (|2Tj) is verified. Then, if we apply k time 
derivatives to the Eq. ([20]) , we straightforwardly obtain 

gk+2 



dt k+2 



Pint(r,t) - p' int (r,t) 



N ( N 



t=t m\N — 



■j") V r- J d^l...dCN-2\^int(^l, V, *o) | 2 Wk (£l , ^ 5 to)) 



(23) 



As, for every physical potential, W^^Wki^i, £jv— lj to) is a real and continuous function in the whole position 
space, and as |^>i n t(£i, £jv-i; *o ) 1 2 is a real and positive function in the whole position space, we can apply the 
integral mean value theorem generalized to many variables functions (demonstrated in Appendix [C| to the previous 
expression. We obtain 

Qk+2 



3(#L, ...,/?at_ 2 ) : m 



Q t k+2 

= V r 
= V r 



Pint(*,t) - p' int {r,t) 



V_N£_.Wfc(/3i, ...,/3 N -2, 
V ~^Wfc(/3i, ...,/3jv_2, 



A r 



iV - 1 

N 
N - 1 



r;to)JV 



- 1 



<2£l...d£/V-2|^mt(£l, • 



jV 



AT- 1 



r;io)f 



r ; to)pi nt (r, to) 



(24) 



To prove the one-to-one correspondence V[uj n t](£i, £jv-i; Pmt(r, t) it remains to show that (|24p cannot vanish 
for w in t and w- nt that lead to the relation (|2"Tj) . Then the internal densities pi nt {v,t) and p' int (r,t) would become 
different infinitesimally later than to- We use the reductio ad absurdum method, in the spirit of Refs. 0, HJ- We 
suppose that (|24[) vanishes, which implies: 



N - 1 
N 

= - f dr 



N \ N 

drw k (Pi, —,0n-2, "77 — rr;t )V r . V jv^iu fc (/3i, ...,/3 N - 2 , — — -r; t )p mt (r, t ) 



Vnt Wk{Pl, ...,/3jv-2, 



A" - 1 
jV 



iV - 1 



AT- 1 



-r;ir 



pint (r, t ). 



(25) 



As Wfe is a many-body function, the Eq. (|2"Tj) does not imply that V(/3i, ...,/3n-2) : V^ N _ 1 Wfc(/3i, Pn-2, £jv-i! *o) 7^ 
in the general case. However, we check if this relation holds for the particular form (H|) we choose for V. 
Inserting the results of Appendix [Al in Q and (|22[) . we obtain, if A/" > 2 (the case A~ = 2 will be discussed later on) 

w/s(/3i, ...,/3jv-2,£iv-i;to) = (26) 



d fc . ,N-1, 



N-2 



Z ^fe Aw ™t(7,-^v-i;*) 



t=t <9i fe ' 



JV-2 



1=1 



where we defined 



"/N-2 



N-2 
N - 1 







JV-2 



JV-2 

and Vi e [1, N - 3] : 7i = - V — — ■ /3 Q 

2 + 1 -^-^ a + 1 



(27) 



a=i+l 



The form of the third term of the right hand side of Eq. ([21))) comes from the fact that X^=i( r « ~ — 0> which 
implies, using the Appendix [Al that — J2a=i STT^ a = — Y^i=i 7»- We see fr° m Eq. ([2"T[) that the set (71, ...,7^-2) 
is perfectly defined by the set (J3i, (3 N-2) an d vice versa. We now can calculate 

^tN-i w k(Pl, ■•■) &N-2, £jV-l5 h) = 



JV-2 



JV-2 



A^ 



A r 



(28) 



where we introduced 



D(r) 



'dt k 



Av int (r;t)\ 



t=t 



(29) 



for simplicity. We now check if 3(@i, pN-2) ■ V^ JV _ 1 iUfe(/3i, Pn-2, £jv-i; *o) = , which is equivalent, according 
to (27]) and (J2SJ), to check if 



JV-2 



JV-2 



3(71, ..,7JV- 2 ) : (AT - 1)D((jV - l)r) = ^ D( 7< - r) + D( - £ 7i - r) 



(30) 



i=i 



i=l 



Some mathematical considerations show that this equa- the set of (71, 7^-2)1 instead if D(r) = const. But 
tion cannot be fulfilled for all r when N > 2, whatever 
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if D(r) = const., then Ai>i„ t (r;t) should be equal to 
b(t).r + c(t)/N, according to (f2"5|) , which is forbidden 
by the condition (|2Tj) . cf. discussion of the i llll Al 

It remains to discuss the case N = 2. It is easy to 
show that then, we have V^ N _ 1 Wk(^N-i't to) = — |Df — 
l&v-i) + |0(|CiV-i), which is null if D(r) is any par 

function of r. But if D(r) is par, then ^ r Awi n t(r; t) 
should be an odd function of r (up to an additional time- 
dependnt function), according to (|29|) . which is also for- 
bidden by the condition ([2"Tj) . cf. discussion of the ^III Al 
Thus, we can conclude that, in our case 

V(/3i, ...,j3 N - 2 ) : V^^WkiPi, ...,/3jv-2,£jv-i;*o) ^ . 

We immediately deduce the incompatibility of this rela- 
tion, which is a consequence of (|2ip and of the particular 
form (|4|) of V, with (|25| . Thus, the hypothesis we made 
is absurd: Eq. (|2^|) cannot vanish if V[A«j n t] ^ c(t), so 
that the internal densities Pi n t(r,t) and p' int (r,t) become 
different infinitcsimally later than t$. As a consequence, 
the map G, defined in (JTHJ) , is invertible (up to an addi- 
tive time-dependent function) and IV'i'ntW) can ^ e wr it- 
ten as a functional of the internal density (we use the 
notation (|17p). Thus, any expectation value of an oper- 
ator O which does not contain a time derivative can be 
written as a functional of pi„t as the phase cancels out. 
This achieves to prove the time-dependent Internal DFT 
theorem (which is a variant of the Runge-Gross theorem 
[1, H[ for self-bound systems and internal densities). 

Mind that all the previous reasonings hold only for a 
fixed initial state ipint(to) (and a given type of particle), 
so that ip® nt is not only a functional of pi nt , but also 
depends on tpi„t(to)- This will be discussed further. 



C. Link with traditional (time-dependent) DFT. 

We stress here the link and differences between the 
traditional DFT and internal DFT potentials. We recall 
that the form of the potential v ex t of traditional DFT can 
be fundamentally justified starting from the laboratory 
Hamiltonian of an isolated molecule where the nuclei are 
treated explicitely. As a molecule is a self-bound sys- 
tem, one can apply the Jacobi coordinates method. We 
denote the N electronic coordinates related to the labo- 
ratory frame as r^, the nuclear cm. coordinate as R™ uci 
and the N electronic coordinates related to the cm. of 
the nuclei as x\ = — R ,mci . A key point concerning the 
molecules is that, as the nuclei are much heavier than the 
electrons, the cm. of the whole molecule coincides with 
R" uci , and it is an excellent approximation to apply the 
Jacobi coordinates to the nuclear coordinates only. As a 
result, the cm. motion will be described by a T(R nucl ) 
wave function. The redundant coordinate problem (thus 
the cm. correlations) will concern the nuclei only, and 
will be "external" to the electronic problem: the N elec- 
trons are still described by N coordinates. Then, if one 



decouples the electronic motion from the nuclear one do- 
ing the clamped nuclei approximation, the interaction of 
the electrons with the nuclear background is described 
by a potential of the form v ext{^i — H nucl ), which 

becomes ^ i « ext (r^) when moving to the cm. frame. We 
then recover the form of the traditional DFT potential. 
The potential v extl which is internal for the (self-bound) 
molecular problem, becomes external for the pure elec- 
tronic problem. Those considerations also hold in the 
time domain, the difference being that the potential 

N 

5> e:ct (r 4 -R"" ci ;t) (31) 

i=l 

can then contain an explicit time dependence in addi- 
tion to the part which describes the interaction of the 
electrons with the nuclear background. We recover the 
traditional time-dependent DFT potential 0, i, HE 
when moving in the cm. frame. 

Those reasonings explicit the link between the tradi- 
tional DFT potential expressed with the laboratory co- 
ordinates, Eq. (|3"Tj). and the Internal DFT potential ex- 
pressed with the laboratory coordinates, Eq. ((T|). They 
both act only on the internal properties, and not on the 
cm. motion (because it is anyway not comparable to 
the experimental one). The difference is that as, in the 
molecular case, some particules are much heavier than 
the other, it is a very good approximation to assimilate 
the cm. of the whole molecule with R ntici , which permits 
to neglect the cm. correlations for the electronic system, 
and to justify the clamped nuclei approximation. This 
simplifies greatly the electronic problem and the tradi- 
tional DFT can be used to study it. When the particules 
constituting the self-bound system have nearly the same 
masses, as it is the case for the nuclei or the He droplets, 
the total cm. (R) should be calculated with all the par- 
ticules, so that the cm. correlations will concern all the 
particules, and no clamped approximation can be justi- 
fied. Then, we should use the formalism proposed here. 

IV. TIME-DEPENDENT INTERNAL 
KOHN-SHAM SCHEME. 

We now provide a practical scheme to calculate the 
internal density pinti which consists in the generaliza- 
tion of the stationary Internal KS scheme of [lj] to the 
time-dependent case. First, we note that for any normal- 
izablc initial state \ipint(to))j which are the only allowed, 
the "internal" Schrodinger equation (fit)]) stems from a 
variational principle on the "internal" quantum action 

Aint = [ dt(ip int (t)\ihdt - H in t\lpint(t)). (32) 

As the function c(t) possibly contained in the potential 
Vi n t is perfectly canceled by the time derivative of the 
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corresponding phase e W *( t )/ B of ipint, see (fl7|) . we have 

A mt = £ dt^LWWI^-^lCtlftntKt)) if Vint 

is chosen so that no additive time-dependent function 
can be split. Thus, the internal quantum action can be 
considered as a functional of pint ■ Its J t 1 dt(tpf nt (t)\ihdt — 

J2a=i JjT U{u]\ip® nt (t)) part is a universal functional 

of pi n t in the sense that, for a given type of particle (a 
given interaction it), the same dependence on pint holds 
for every V[v lnt ], thus v int (see (@}). 

Using the Eq. © , we develop the "internal" quantum 
action as 



Aint[Pint] 



dt(^ nt (t)\^d t 



N-l 

E 

a=l 



where we used to obtain the penultimate equality. 

We see that the potential J2iLi v int(*i — R;i) that is 
N body with respect to the laboratory coordinates (and 
(N — 1) body when expressed with Jacobi coordinates), 
becomes one body (and local) when expressed with the 
cm. frame coordinates (mind that pi n t is defined in the 
cm. frame, i.e. that r is measured in the cm. frame, cf. 

Applying ([34]) to the second term of the action in- 
tegral (33J) gives (^ nt (*)I^M(6, 6r-iM„ t (*)) = 
\ Jdr dr'j int (r, r'; t)u(r — r'), where % nt is defined in 

The action integral (|33|) can thus be rewritten 



2p a 



(33) 



dt(^ nt (t)\UM(^-^N-lM nt (t)) 



A int [ Pint ] = / dt(^ nt mm-Y,2^\^t(t)) 



ti 



dt(^ nt (*) \V[v int ] U'l i\ ,;f) | (*) ) ■ 



'to 

1 

2 



eft drdr'jint{r, r';t)u(r 



To rewrite its last two terms in a more convenient 
way, we establish a useful relation. For any function 
/(ri, ...,r?v;t) of the laboratory coordinates, expressible 



dt dr Vint(r;t)pint{r,t). (36) 



with the Jacobi coordinates [we denote F(£i, 
we have 



«t(*)|f(fc 
#i • • • 



^ JV _ lJ F(a,... ! ejv-i;t)|V'L(6 J 



■dCiv-i<5(R)^(e 



dri • • -drAr(5(R)/(ri, . 



&v-i;i) 

x|^L(Ci,- 

r w;<)kL( r i: 



,6v-i;*)L 



(34) 
> ^jv— i; *) I 



■,^N-i;t)\ 
..,r N ;t)\ 



Up to now we did not do any hypothesis. To recover 
the Internal time-dependent KS scheme, we assume, as 
to obtain the traditional time-dependent KS scheme [J, 
5], that there exists, in the cm. frame, a ./V-body non- 
interacting system (i.e. a local single-particle potential 
vs) 



h 2 A 
2m 



+ v s (r,t))4 nt (r,t)=iW t ipi nt (v,t) (37) 



which reproduces exactly the density p% n t of the interact- 
ing system (mind that p^t is defined in the cm. frame) 



N 



Pint 



We see that the "internal mean values" calculated with 
ipint expressed as a function of the (N— 1) coordinates £ a , 
can also be calculated with ipint expressed as a function 
of the N coordinates r^. As one of them is redundant, a 
<5(R) which represents the dependence of the redundant 
coordinate on the others appears [3fj| . 
The relation (f3"l| leads to 



(38) 



(V> 



(t)\v[v m Mi 

dri 



N 



?«*(*)) 



dr N S(R)^v in t(ri 



R;0lC f (ri, 



•,rjv;t)| 



; = 1 



N 

E 

i=i 



N 



dr v int (r;t) / dr x 



i=l 



^2 dr v int {r;t) 



x\ip? nt (r u ... 
Pint{r,t) 



dr N 6(R) 

,T N ;t)\ 2 8(r-{Ti-R)) 



N 



dr v int (r;t) p mt (r,t), 



Even if only (TV — 1) coordinates are sufficient to describe 
the internal properties, they still describe a system of ./V 
particles. Thus, we have to introduce ./V orbitals in the 
KS scheme (as we did) if we want them to be interpreted 
(to first order only) as single-particle orbitals and obtain 
a scheme comparable (but not equivalent) to mean-field 
like calculations with effective interactions. 

In (|37[) we implicitely supposed that the particles are 
2 Fermions (a KS scheme to describe Boson condensates 
can be set similarly equalling all the f\ nt ). Unique- 
ness of the potential i>s(r, t) for a given density Pmt(r, t) 
(and initial |v?i„ t (io)) which yield the correct initial den- 
sity pint (r, to)) is ensured by a direct application of the 
traditional time-dependent DFT formalism [3, [f|- Of 
course, the question of the validity of the KS hypothesis, 
known as the non- interacting v-representability problem, 
remains, as in traditional (time-dependent) DFT [l], Q. 

To use similar kinds of notations than the traditional 
DFT ones, we add and substract to the internal action 
^ 'integral (f3"6")) the internal Hartree term 
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A H [pint] = \\ t l dt Jdrdr' p int (r, t) p int (r' ,t)u(r - r'), 
the non-interacting kinetic energy term 

ntvOl^lvLW) and tllC 

int (t)\ihd t \ipl nt (t)) term. This permits to 
rewrite the "internal" action integral (f36|) as 



-A X c[pint] 



dt / dr w mt (r; i) p int (r, i) (39) 



<E 
t=i 



2m 



|^ rat (i)) — Ai3-[pi„t] where the internal exchange-correlation part is defined as 



Axc[pi 



- [ dt [drdr'(^i nt (r,r';t)-pi n t(r,t)pint(r',tf)u(r-r') 



1 N-l 2 N 2 



A' 



i=l 

I 



(40) 



We see that it contains the exchange-correlation which 
comes from the interaction u (first line of (|4"01) ). but also 
the correlations contained in the interacting kinetic en- 
ergy (second line of (00])) and in the interacting "ihd t " 
term (third line of (|40|) ). A key point is that, as the KS 
assumption implies ipl nt [p int ] [3,11, H], A X c[Pint](t) can 
be written as a functional of pi„t (for given \ipi nt (to)) 
and {|</?!„t(£o))} which yield the same initial density 
Pint(r,t )). 

It remains to vary the "internal" quantum action (|39|) 
to obtain the equations of motion (which define pint)- 
Vignale, see Ref. [32j |. showed recently that the correct 
formulation of the variational principle is not to station- 
arize the quantum action, i.e. 5Ai n t[pint] = as done so 
far 0,1, EH, but 

SA int [p in t} = ityint\Pint](tl)\6i>int\Pint](tl)) 



> (i>lt [Pint](tl ) 1 8^ nt [ Pint ](ti )) (41) 



(where tpf nt is the Slater determinant constructed from 
the tpl nt ) ■ The two formulations lead to identical final re- 
sults for theorems derived form symmetries of the action 
functional because compensations occur [32| . but Vig- 
nales's formulation permits to solve the causality paradox 
of the previous formulation. 

Varying ([4Tj) with respect to the <p-* t (r, i), with t 6 
[t0)*i]i leads straightforwardly to the Internal time- 
dependent KS equations for the (p\ nt 



I % \ ■ 

(42) 

with the potentials 



U H [pint]{r,t) = 
Uxc[Pint}{r,t) -- 



SAh [pint] 

5p int (r,t) 

_ SA XC [pint] 

Sp int {r,t) 



i faint [Pint]{h) 



Hint [Pint]{h) ~ 

5p int (r,t) ' 



+ i{^int[Pint]{h) 



Spi n t(r,t) 



) (43) 



which are local as expected (vs = Uh [pint] + Uxc[pint] + 
v^t with the notations of Eq. (1371) ). Note that the varia- 
tional formulation of Vignale [32| leads to the addition of 
the last two terms in the definition of Uxc[pint](r,t), see 
Eq. (|4"3"| . compared to the traditional result obtained by 
stationarization of the action. It are those terms which 



permit to solve the causality paradox [32j. 

Equations (|4"2"| have the same form as the tradi- 
tional time-dependent KS equations formulated for non- 
translationally invariant Hamiltonians [3, 1, 0] and per- 
mit to define pi„t through (|38[) . Here, we have justified 
their form in the cm. frame for self-bound systems de- 
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scribed with translationally invariant Hamiltonians. one can show that the interacting kinetic energy term 

But there is a major difference with the traditional and the interacting "ihd t v term can be rewritten [30| 
DFT formalism. Following similar steps as in Eq. (|34|) . 



2V-1 9 n 



a— 1 2—1 

(Ct(uKW°nt(u) = / dr 1 ..-dr JV *(R)Ct(ri»»-»r J v;WCt('i»»-»rjv;*), (44) 



which makes it clear that the differences 
with the non-interacting kinetic energy term 
J2 l= iJ dripi* nt {r)l^<pl nt {r) and the non-interacting 

SM) t term" E^i(Vint(*)|i^lvU(*)) ( found in thc 
exchange-correlation functional (|40|) ) come, on the one 
hand, from the correlations neglected in the traditional 
independent-particle framework, but also from the cm. 
correlations described by the S(R) term in (|44|) . which 
does not appear in traditional time-dependent DFT 
[J, [f|. The inclusion of thc cm. correlations in the 
exchange-correlation functional (|40[) and potential (|43)) 
is the main difference with the traditional KS scheme, 
and is a key issue for self bound-systems as atomic 
nuclei. 



Mind that all the previous considerations only hold 
for fixed initial states \ipint(to)) and {IVint^o))} which 
should of course give the same initial density /Oj n t(r, to) 
(and also for a fixed type of particle). As a consequence, 
■0° nt is not only a functional of Pi„t, but also depends 
on the initial state \ipi n t(to)), an d Uxc, Eq. also 
depends on the initial orbitals {I'/'int^o))}- An impor- 
tant difference to the ground state Internal DFT formal- 
ism / KS scheme presented in [l4| is that |Vw(*o)) an d 
the {\<-p\ nt {to))} cannot necessarily be written as func- 
tionals of pi„ t (r,t ). However, as underlined in 0, Q, if 
one starts from initial states \ipint(to)) and {\fl nt (to))} 
that are non-degenerate ground states, i.e. that can be 
written as functionals of pi n t(r,to) [14| . tpi nt and Uxc 
become functionals of pi n t(j,t) alone. Then, in the limit 
of stationary ground states, the theory reduces to the 
stationary Internal DFT / KS. 



We recall that, as in traditional DFT, the previously 
discussed functionals are defined only for internal den- 
sities p^t which correspond to some internalpotential 
Vint , called v-representable internal densities [4|, [5|]. Up 
to now, we do not know exactly how large the set of v- 
representable densities is. This has to be kept in mind 
when variations with arbitrary densities are done, as to 
obtain the time-dependent KS equations. 



V. CONCLUSION. 

In summary, we have shown that, for a fixed initial 
state, the internal wave function, which describes the in- 
ternal properties of a time-dependent self-bound system, 
can be written (up to a trivial phase) as a functional 
of the internal density. This implies that the "inter- 
nal" expectation values of any observable (which does 
not contain a time derivative), that are of experimental 
interest, can be regarded as functionals of the internal 
density. Then, we set up, in the cm. frame, a practi- 
cal scheme which permits to calculate the internal den- 
sity and whose form is similar to the traditional time- 
dependent KS equations, the difference being that the 
exchange-correlation functional contains the cm. corre- 
lations. 

This work is a first step towards the justification to 
the use of density functionals for time-dependent nu- 
clear mean-field like calculations with effective interac- 
tions [l^, [24[ , proving that there exists an ultimate func- 
tional which permits to reproduce the exact internal den- 
sity (up to the non-interacting v-representability ques- 
tion). If this functional was known, there would be no 
need for a cm. correction. 

Practically speaking, the time-dependent Internal KS 
scheme can describe, for instance in the nuclear case, 
the collision of two nuclei in the frame attached to the 
total cm. of the nuclei. Then, Vint is zero but the de- 
pendency to the initial state allows to start from a state 
which corresponds to two nuclei with different velocities, 
or "boosts" (choosen such as the total kinetic momentum 
is zero because we are in the cm. frame). According to 
the choice of the boosts, we can describe a wide variety of 
physical phenomena, from nuclear fusion p3j to Coulomb 
excitation (33| . One of the nuclei can also simply consist 
in a particule as a proton, to describe the excitation of a 
nucleus by diffusion. 

A case where a non-zero Vi n t would be interesting could 
be the case of the laser irradiation (vi n t would then con- 
tain a laser potential switched on at t > to). This is not 
of major interest in the nuclear case because, experimen- 
tally speaking, we do not yet have lasers that are suited 
to the study of the laser irradiation of a nucleus. How- 
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ever, this could be interesting in view of a generalization 
of this work to the whole molecule (following from the 
generalization to different types of particulcs, which is 
underway) . 

Many questions remain open. In particular, the ques- 
tion of the form of the potential which describes the 
cm. correlations; in addition to its practical interest, 
this question would also give interesting arguments con- 
cerning the non-interacting v-reprcscntability question. 
Generalization to different types of particles (Fermions 
or Bosons) appears desirable. Finally, the same reason- 
ing should be applied to rotational invariance to formu- 
late the theory in term of the so-called "intrinsic" one- 
body density [3J| (which is not directly observable) . This 
is more complicated because rotation does not decouple 
from internal motion, but it should be interesting con- 
cerning the symmetry breaking question. 



Vi € [2, N - 2] : 

• _ 1 N-2 

r i - R = - — rCa - -776V-I, 

l a + 1 N 



r l — R = - ^ — £ a - — fjV-1. 

^— ' a + 1 N 

a—l 

The £n-i is the only Jacobi coordinate to appear in the 
expressions (in function of the Jacobi coordinates) of all 
the {r^ — R}. Thus, all the Vi n t(ri — R;£) terms that 
enter into the calculation of the total internal potential, 
Eq. (j4|), will contain £jv— i- 



Acknowledgments 

The author is particularly grateful to M. Bender, 
E.K.U. Gross, and E. Suraud for enlightening discussions 
and reading of the manuscript, and thanks the referee for 
his pertinent remarks. The author thanks the Centre 
d'Etudes Nucleaires de Bordeaux-Gradignan for warm 
hospitality, and the Institut Universitaire de France and 
the Agence Nationale de la Recherche (ANR-06-BLAN- 
0319-02) for financial support. 



APPENDIX A: EXPRESSION OF THE {r^ R} AS 
FUNCTIONS OF THE {£ Q } COORDINATES. 



APPENDIX B: EXPRESSION OF THE 
INTERNAL PROBABILITY CURRENT j ini (r,t). 

The one-body total laboratory density and probability 
current are defined as 

p(r,t) = N J dr 1 ...dr N - 1 \ip(r 1 , r w _i, r; t)| 2 

j( r >*) = n — :N / dr 1 ...dr N ^ 1 ip*(r 1 , rjy-i, r; i) 
Zmt J 

x V r ^(ri,..., rjv-i, r; t) + c.c. (Bl) 



Using the relations J3]), one could show that the {r.; — 
R} can all be written in function of the coordinates. 
We obtain, as a result 

N-l 
rjv — R = — — — sat-1) 

N-2. 1 
r w _! — R — ^37^-2 - jy6v-i, 



R 



N-3 



N 



1 



N - 1 



They satisfy the "laboratory" continuity equation 



5 t p(r,t) + V r .j(r,t) =0. 



(B2) 



We do some manupulations on the laboratory probabil- 
ity current j using the Jacobi coordinates (c.c. denotes 
the complex conjugate; the second equality is obtained 
using V rN = V^ N _ X + Vr/7V, by definition of the Jacobi 
coordinates) 



j(r,t) = -j—N jdr 1 ..dr N -idr N 5(r-r N )^*(r 1 ,..,,r N _i,r N it)V rN Tp(r 1 ,.,.,r N -i,r N ]t)+c.c. 

= ^/^■■•^ 1 (iv^) 3 <^- 1 -F^ (r - R) ) 

x r* (R, t)^* nt (f 1 , ■ ■■ , 6v-i ; t) (v eN _, + r(R, t)ip int (&,..., £w-i; t) + c.c. 

= J dR|r(R,t)| 2 x _tv(__) y da---^jv-2^„ t (a,---,^-2^;OV,^(a,---,^-2^;t)[ = ^ (r 
+ i / dRr *( R ^)^r(R,t) xN(-^y J dei...de^-a|^t(6,-.fjv-a )] ^Y(r-R);t)| 2 + c.c. 
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J dR|r(R,t)| 2 j mt (r-R, t)+ J dRp int (r-R, t)j r (R,t), 



(B3) 



where we introduced the internal one-body density pip , the cm. probability current jr(R, t) = 
2jyjT*(R, £)VrT(R, t) + c.c. (M = Nm is the total mass), and the internal probability current ([13]) . The mean- 
ing of Eq. (|B3|) is clear: the laboratory probability current is the sum of the cm. probability current and of the 
internal probability current, both convolued respectively with the internal one-body density and the cm. one-body 
density. One can show that jr and jj n j satisfy both independent continuity relations. It is trivial, using (J5]) and ([5]), 
for jr. For ji„i, we calculate (with the help of ([6]) and (|T0|) ) 



N 



N - 1 



d t I d£i...d£/v-i£tev-i 



N 



N - 1 



1 f N \ 3 f N h 2 

-fr\ N _ I ) I dCi-d^N-iS^N-! - — — T r)ipi nt (£i, ...^N-^t)-^— — A ilN _ 1 ip m t(ti, ...,£jv-i;i) + c.c. 



N- 1 



2/ijV- 



h 

.JL( 



N 



2mi \N - 1 



— ( N ~ 1 ) j d^i..^ N - 2 4>tnt(£li---^N~2,V]t)A y -4; m t((,l,...,£N-2 1 V- 1 t) 

( N \ 3 /" 

_ 1 ) V ^ / d £l---dtN-2ll>tnMli -! &V-2, ^5 t)Vvtpi„t{£l, -,£n-2, V] t) 

V r . / dCl...<^iV-2V , *nt(6, ■•• ! ^JV-2,^*)V„'!/'int(Cl) &V-2, ^5 *) 

I 



+ C.C. 



+ c.c. 



+ c.c. 



(To obtain the last equality we used the fact that, by def- 
inition, [Xn-i = ^w^ m -) From this relation we deduce, 
using (|lip and (p~3|) . the "internal" continuity equation 

d t p mt (r,t) + V r .ji„ t (r,i) = 0. 

This reinforces the interpretation of jmt as the internal 
probability current. 

APPENDIX C: INTEGRAL MEAN VALUE 
THEOREM FOR FUNCTIONS OF MANY 
VARIABLES. 

We give the generalization of the mean value theorem 
[35| to functions of an arbitrary number of variables. One 



starts from two real functions of A variables 

V(a?i, ...,xa) G 5fte : / : [x x , ...,x A ] i-» 

g : [xi, ...,xa] i-> 3?e. 

We suppose that they are integrable in a domain D, that 
/ > in D, and that g is continuous in D. We define 

m = inf{g{xi 7 ...,x A ); (x%, ...,x A ) G D} 
M = sup{g{x 1 ,...,x A )\ (xi, ...,x A ) G D}. 

As / > 0, m/ < /g < Mf, which we integrate 



ml dxi...dx A f(xi, ...,x A ) < I dx\...dx A f{xi, x A )g(xi, x A ) < M dxi...dx A f(x 1 , ...,x A ) 

JD JD JD 

3C £ [m, M] : / dxi...dxAf(xi, x A )g(xi, xa) = C I dx 1 ...dx A f(x 1 ,...,x A ). 

J D J D 

As we supposed that g is continuous in D, we deduce that 3(x[, ■■■,x' A ) G D : g{x' 1 , ■■■,x' A ) = C, which implies 
3(x[, ...,x' A ) G D : I dx 1 ...dx A f(x 1 ,...,x A )g(x 1 ,...x A )=g(x' 1 ,...,x' A ) dx 1 ...dx A f(x 1 ,...,x A ). 

JD J D 
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